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Abstract 

Emission from the 2i2 - In line of H2GO has been detected and marginally resolved toward LkGa 15 
by the Nobeyama Millimeter Array. The column density of H2GO is higher than that observed in DM Tau 
and than predicted by theoretical models of disk chemistry; also, the line-intensity profile is less centrally 
peaked than that for GO. A similar behavior is observed in other organic gaseous molecules in the LkCa 
15 disk. 

Key words: ISM: molecules — stars: individual (LkCa 15) — stars: pre-main-sequence — circumstel- 
lar matter 


1. Introduction 

It is now well-established that more than 50% of young 
T Tauri stars are surrounded by a disk of circumstellar 
material. These disks are important because they are the 
birth sites of planetary systems. The spectral energy dis¬ 
tributions indicate that disk masses are in the range of 
~ 10“^ - 1O“^M0, which is consistent with the solar neb¬ 
ula model proposed for the origin of our own solar sys¬ 
tem (e.g., Beckwith, Sargent 1996). Interferometer ob¬ 
servations in the dust continuum and CO line emission 
have spatially resolved some of these disks (e.g., Saito et 
al. 1995; Dutrey et al. 1996; Hogerheijde et al. 1997; 
Guilloteau, Dutrey 1998; Mundy et al. 2000; Duvert et 
al. 2000), which are found to have radii of ~ 100 - 800 
AU. However, our quantitative understanding of the phys¬ 
ical properties of such disks, e.g. their radial and vertical 
temperature and density structures, or their gas survival 
timescales, is still poor. Another important question is 
the chemical evolution of gas and dust as it is transported 
from the interstellar medium to the interiors of circumstel¬ 
lar disks, and the impact of this chemistry on the nature of 
icy planetesimals, such as comets and Kuiper Belt objects 
(van Dishoeck, Blake 1998; Ehrenfreund et al. 1997). 

We are involved in a project to systematically investi¬ 
gate the physical and chemical properties of circumstellar 
disks. Single-dish spectra of simple molecules such as CO, 
HCO'*', CN, HCN, CS, and H2CO have been obtained at 
the IRAM 30 m telescope, JCMT and CSO (van Zadelhoff 
et al. 2001; Thi 2002, Thi et al. in preparation). The ra¬ 
tios of high-J and low-J lines indicate that the detected 


gaseous molecules reside in a low-temperature (e.g. ~ 30 
K) region, which corresponds to disk radii beyond ^ 100 
AU, assuming the temperature distribution in typical disk 
models. In addition to single-dish observations, interfer¬ 
ometer observations of low-J transition lines of several 
molecular species have been performed by Duvert et al. 
(2000) and Qi (2001). Although interferometer observa¬ 
tions of molecular lines in disks are still rare because of the 
faintness of the objects, they provide an important probe 
of the physical and chemical gradients in the disks; the 
line intensities depend on the density, temperature, and 
molecular abundance, which can be affected by various 
chemical processes in each region of the disk. 

In this paper we report on observations of the H2CO 
2 i 2 - 111 line at 140.84 GHz toward LkCa 15 using the 
Nobeyama Millimeter Array (NMA). LkCa 15 is a rela¬ 
tively old (~ 1 X 10^ yr) solar-mass T Tauri star located 
in the outer regions of the Taurus molecular cloud at ~ 140 
pc. Its disk mass is estimated to be ~ 0.02 - O.O6M0 from 
dust continuum observations (i.e. Thi 2002; Kitamura et 
al. 2002). The formaldehyde molecule is of special inter¬ 
est, because its line ratios are sensitive to the temperature 
and density (e.g. Mangum, Wootten 1993). Single-dish 
observations of various H2CO lines have been obtained by 
Thi (2002) using the IRAM 30 m telescope and JCMT. 
It is also one of the most complicated molecules detected 
so far in disks, serves as a probe of grain-surface chem¬ 
istry, and can rapidly polymerize to create complex or¬ 
ganic solids under appropriate conditions (Schutte et al. 
1993). Observations of the transitions in the 1.3 mm at¬ 
mospheric window using OVRO and a comparison with 
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the data presented here will be reported in a forthcoming 
paper. 

2. Observations 

We observed the 2i2 ~ In line of H 2 CO toward LkCa 
15 with the NMA in 2000 - 2001. The results presented 
here are based primarily on 16 hr of array measurements 
using six antennas in the NMA low-resolution (D) config¬ 
uration, under clear-sky conditions. Similar quality data 
were obtained over a 6-hr period in the high-resolution 
(AB) configuration, but no line emission was detected. 
The beam sizes of the D and AB configurations are about 
4."5 and l."2, respectively. 

A dust continuum emission at A = 2 mm was detected 
with a good signal-to-noise ratio with the NMA SIS double 
sideband (DSB) receivers. The zenith DSB system tem¬ 
peratures during the observations were typically 200 K. 
For back ends, the Ultra Wide Band Correlator (UWBC, 
Okumura et al. 2000) and a high-dispersion FX cor¬ 
relator were operated simultaneously. Phase-switching 
techniques were used to separate the continuum visibil¬ 
ity data for the lower (128.840 ± 0.512 GHz) and up¬ 
per (140.840 ± 0.512 GHz) sidebands from the UWBC. 
Spectral line visibility data collected by the FX correlator 
were obtained by subtracting the continuum level, which 
was estimated by averaging the line-free channels. 

The response across the observed passband for each 
sideband was determined from 40-min observations of 3C 
454.3. The gain calibrator 0446-1-112, used to determine 
the flux density scale, was observed every 20 min. The 
2 mm flux density of 0446-1-112 was estimated to be 1.57 
Jy at the time of these observations from a comparison 
with the visibilities measured for Uranus (Griffin, Orton 
1993), with an overall uncertainty in the flux calibration 
of ^^10%. Final images were made using the calibrated 
visibilities within the AIPS package developed at NRAO. 

3. Results 

3.1. Detection of H 2 CO in the Disk of LkCa 15 

The continuum emission from the LkCa 15 disk was 
clearly detected; an image using only the low-resolution 
configuration upper sideband data is shown in figure 1. 
The total dust flux at 141 GHz is 25 mJy and is unre¬ 
solved, leading to an upper limit to the dust disk radius 
of ~ 550 AU. 

The line spectrum for a 10" box centered on the stel¬ 
lar position is shown in figure 2. Emission was detected 
in the range ±2.5 km s“^ around the stellar velocity 
Vs = fLSR = 6 km s“^. Since the signal-to-noise ratio is 
not high enough to draw a image with the velocity bin of 
~ 0.7 km s“^ adopted in figure 2, we integrated the data 
from u — Us = —2 to ±2.6 km s“^ to obtain figure 3. The 
contour interval is 21 mJy beam“^, which corresponds to 
the 1 a r.m.s. noise level. The position of the H 2 GO 
emission coincides with that of the dust continuum. It 
should also be noted that previous CO observations show 
no foreground or background molecular gas towards the 
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object (Duvert et al. 2000). We therefore conclude that 
the H 2 CO emission arises from the disk around LkCa 15. 
The line spectrum is consistent with the double-peaked 
shape characteristic of Keplerian disk velocity fields, al¬ 
though we could not detect any gradients in the channel 
maps with higher kinematic resolution (e.g. 1 km s“^) 
because of the low signal-to-noise ratio. 

The H 2 CO disk is marginally resolved; it is more ex¬ 
tended than the continuum image, which is typical for 
molecular line images because of the increased opacity 
as compared to that of the dust (Dutrey et al. 1996). 
The H 2 CO disk elongation from north-west to south-east 
is caused by the disk inclination from the line of sight. 
Interferometer maps of CN line emission show similar 
elongation along the same direction, and CO observations 
reveal a velocity gradient that runs along this major axis 
(Duvert et al. 2000; Qi 2001). The radius of the H 2 CO 
disk, measured at the 2-cr level, is about 650 AU, which 
is similar to the CO disk radius deduced from Plateau de 
Bure interferometer observations (Duvert et al. 2000). 

The total flux density of the H 2 CO line averaged over 
the 4.6 km s“^ interval at the stellar velocity is about 180 
mJy, leading to an integrated line intensity of 820 mJy km 
s“^. This is similar to the value obtained from single-dish 
data at the IRAM 30 m telescope by Thi (2002), which 
is about 664 mJy km s“^, as is the estimated line width 
of ~ 4 km s“^. The interferometer observations therefore 
trace the same gaseous component detected by the IRAM 
30 m telescope. 

3.2. Molecular Column Density 

The peak flux of the H 2 CO line is about 95 mJy 
beam“^, which corresponds to an antenna temperature 
{Tpf) of 0.28 K for a 4."48 x 4."13 synthesized beam. The 
averaged antenna temperature of the line within the 2 a 
contour is about 0.16 K. Considering typical disk models, 
because the gas density is higher than the critical den¬ 
sity of the line, ~ 10^ cm“^, the line is almost in local 
thermodynamic equilibrium (LTE). Since the disk image 
is marginally resolved, it is not likely that the beam-filling 
factor is much smaller than unity. If the line is optically 
thick, with the above two conditions, the antenna tem¬ 
perature should be close to the kinetic temperature of 
the molecular gas. The observed antenna temperature, 
however, is much lower than the kinetic temperature of 
molecular gases in typical disk models, ^ 20 K. Hence, 
the molecular line is likely to be optically thin. In this 
subsection we derive the molecular column density in the 
disk using LTE models. 

The averaged antenna temperature, 0.16 K, leads to 
an averaged upper level column density of ~ 5.4 x 10^^ 
cm“^ from the formula AAip = 87rfc^^TaAu/(/ic^Aui), in 
which Aui = 5.3 x 10“® s“^ is the Einstein A coeffi¬ 
cient (e.g. Goldsmith, Langer 1999). The total col¬ 
umn density of ortho H 2 GO can be calculated from 
Ntot = A^upQ/[ 5 up exp(—U/fcT)], where Q is the rota¬ 
tional partition function, = 15 is the degeneracy of 
the upper state, and E = 21.9 K is the upper state en¬ 
ergy. Theoretical models of disk chemistry predict gaseous 
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Fig. 1. Integrated intensity dust continuum map of the cen¬ 
tral 60'^ region towards LkCa 15. The contour intervals are 
2 (T, starting at ±2 cr (1 cr = 1.173 mjy beam~^). Negative 
contours are shown as dashed lines. 


Fig. 3. Intensity map of H 2 CO 2i2 - In integrated from 
V — Vs = —2 to -1-2.6 km s~^. The contour intervals are 1 
cr, starting at dil cr (1 cr = 21 mJy beam“^). Negative con¬ 
tours are shown as dashed lines. The cross indicates the peak 
position of the dust continuum. 



v-Vs [km s-i] 


Fig. 2. Line spectrum of H 2 CO 2i2 - In in the central 10^' 
region of the map. The stellar velocity Vs is 6 km s“^. 

organic molecules exist predominantly in regions where 
T ^ 20 K (section 4), while recent disk radiative trans¬ 
fer models predict that the surface region of the disk can 
be heated to temperatures of ~ 50 K by the central star 
(Chiang, Goldreich 1997; D’Alessio et al. 1998). For Trot 
values of 20 or 50 K, the column density of ortho H 2 CO is 

5.4 X 10^^ cm“^ or 1.1 x 10^^ cm“^, respectively, in LTE. 
The unknown ortho-para ratio is another uncertainty in 
estimating the total column density of H 2 CO. In interstel¬ 
lar clouds the formaldehyde ortho-para ratio ranges from 

1.5 - 3.0 (e.g. Mangum, Wootten 1993). Hence, the total 
H 2 CO column density averaged within the 2 a contour, 
which corresponds to a region of about 970 AU x 1250 
AU, is estimated to be 7.2 x 10^^ — 1.9 x 10^^ cm“^. 


4. Discussion 

DM Tail is another T Tauri star/disk system toward 
which H 2 CO has been detected (Dutrey et al. 1997). The 
estimated molecular column density is ~ 1 X 10^^ cm 
which is much smaller than that in LkCa 15. Observations 
of other molecular lines also show that the column den¬ 
sities of gaseous organic molecules are higher in LkCa 15 
than in DM Tau by a factor of ~ 10 (Qi 2001). 

The molecular abundances in the outer regions 100 
AU) of disks have been theoretically investigated by sev¬ 
eral groups (Aikawa, Herbst 1999, 2001; Willacy, Langer 
2000; Aikawa et al. 2002; van Zadelhoff et al. 2003). In 
all such models, the chemistry is intimately tied to the 
physical state of the gas and the properties of the disk. 
For example, the upper panel of figure 4 shows the verti¬ 
cal and radial distribution of H 2 CO in the model by van 
Zadelhoff et al. (2003). They adopted the disk model of 
D’Alessio et al. (1999): solved the ultra-violet (UV) 

radiation transfer and chemical reaction network, which 
includes gas-phase chemistry, freeze-out and evaporation, 
but no active surface chemistry. Near the disk surface, 
molecules are dissociated by UV radiation from the inter¬ 
stellar field and central star. In the disk midplane organic 
molecules are heavily depleted onto grains, since the den¬ 
sities are very high and the temperature is lower than 20 
K, the sublimation temperature of the dominant carbon 
reservoir, CO. Thus, in the outer reaches of circumstel- 
lar disks, gaseous molecules exist in high abundance only 
in an intermediate height layer that is heated by the in¬ 
frared radiation from the surface layer yet shielded from 
the harsh UV radiation field. Specifically, formaldehyde 
is found to have a maximum gas-phase abundance in gas 
at Av ~ 0.4 - 10 mag from the disk surface. 
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Model comparisons with observations require knowl¬ 
edge of the molecular column density versus distance from 
the central star, and the lower panel of figure 4 shows the 
H 2 CO column density calculated by van Zadelhoff et al. 
(2003). Their result is in reasonable agreement with ob¬ 
servations of DM Tau, but smaller than that needed to re¬ 
produce the LkCa 15 observations. Clearly, there must be 
some difference in the physical structure or evolutionary 
paths of DM Tau and LkCa 15 which causes the chemi¬ 
cal variation seen in these two objects. Interestingly, the 
models predict that the molecular column densities do not 
significantly depend on either the disk mass or age; that 
the mass of the intermediate molecular layer does not de¬ 
pend on the total disk mass as long as the disk is as thick 
as Av ^ 20 mag, and that the chemical timescale in the 
intermediate layer is relatively short, ~ 10“^ - 10® yr. On 
the other hand, X-rays from the central star may cause 
variations in the molecular abundances, since the X-ray- 
induced UV destruction of CO can supply the additional 
carbon needed to form other species. Neither DM Tau 
nor LkCa 15 has been detected by ROSAT X-ray sur¬ 
veys. However, considering the high time variability of 
X-rays from T Tauri stars, reobservation of these objects 
with Chandra would be highly desirable. Dust coagula¬ 
tion and sedimentation could also cause variations of the 
molecular abundances since such processes would signifi¬ 
cantly change the temperature distribution and UV radi¬ 
ation field within the disks. Indeed, recent studies on the 
spectral energy distribution indicate the coagulation and 
settling of dust in several T Tauri and Herbig Ae/Be disks, 
including LkCa 15 (Chiang et al. 2001; D’Alessio et al. 
2001). Finally, active grain surface chemistry with evapo¬ 
ration may enhance the H 2 CO abundance, and may cause 
variations of the molecular abundances among objects, if 
they have different temperature or grain size. 

The spatial distribution of molecules is another interest¬ 
ing issue concerning the chemical and physical structure of 
protoplanetary disks. Theoretical models predict that the 
H 2 CO column density is fairly insensitive to the total H 2 
column density, and thus to the disk radius. On the other 
hand, CO is one of the most volatile species after molecu¬ 
lar hydrogen, with a sublimation temperature of only ~ 20 
K, which is reached in the midplane at i? ~ 100 AU. This 
causes a steep gradient in the radial distribution of the 
CO column density. In our observations, the H 2 CO line 
intensity is less centrally peaked than that of CO, which 
was observed by Duvert et al. (2000). However, we can¬ 
not exclude the possibility that the flat distribution of the 
H 2 CO emission intensity is caused by the low spatial reso¬ 
lution and/or low signal-to-noise ratio of our observations; 
the central region of the Keplerian disk corresponds to the 
wing component of the line spectrum, which is more eas¬ 
ily affected by noise than velocities corresponding to the 
outer regions of the disk. 

In order to test these possibilities, we observed the LkCa 
15 disk with the high-resolution (l."2) configuration of the 
NMA for 4 nights. Unfortunately, clear-sky conditions 
were present for only 6 hr, leading to a noise level of 82 
mjy beam“^ or 3.1 K. Since the peak antenna tempera- 



Fig. 4. Distributions of abundance and column density of 
H 2 CO predicted by van Zadelhoff et al. (2003). It is their 
model with Spectrum B; including excess UV radiation from 
the central star with a spectrum similar to that observed 
for TW Hya. The gray scale in the upper panel shows re¬ 
gions in which H 2 CO abundance relative to hydrogen nu¬ 
clei is < 5 X 10“^^ (white), 5 X 10“^^ — 1 X 10“^^ (the dark¬ 
est gray), 1 X 10“^^ — 1 X 10“^^, 1 X 10“^^ — 1 X 10“^° and 
1 X — 2.5 X (the lightest gray). The line contours 

depict the temperature distribution. The lower panel presents 
the radial distribution of the H 2 CO column density integrated 
in the vertical direction. 

ture of 0.28 K in the D configuration corresponds to 3.6 
K in a l."2 beam, these high-resolution observations do 
not add any significant constraints on the spatial distri¬ 
bution of H 2 CO. Some quantitative constraints are given 
as follows. If the temperature of the H 2 CO at i? = 150 
AU is similar to that at i? = 600 AU, the ratio of the 
total H 2 CO column density at these radii should be pro¬ 
portional to that of brightness (i.e. mean antenna tem¬ 
perature), and is less than 3.1/0.16 ~ 19. If the molecular 
layer in the inner radius has a higher temperature than 
that at the outer radius, we obtain a weaker constraint. 
Assuming temperatures of 50 K and 20 K at i? = 150 AU 
and 600 AU, respectively, the column density ratio is less 
than 40. Observations with higher sensitivity and higher 
spatial resolution are therefore critical for an improved 
understanding of the molecular distributions and chemi- 
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cal processes within protoplanetary disks. 
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